Maximal qubit tomography* 
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A new qubit tomography protocol is introduced, based on a continuous positive operator valued 
measure, which is supported by the set of pure states, and equivariant under the symmetry group 
SO(3, R) of the qubit state space. Thus the sample data for the tomography protocol is a discrete set 
of points on the Bloch sphere obtained by a series of independent measurements on identical replicas 
of the qubit. Although the maximum likelihood estimator uniquely exists, it cannot be explicitly 
calculated. Instead, we present another convenient, simple estimator, show that it is unbiased and 
calculate its variance. 

Furthermore, a detailed account is given about known properties of related discrete qubit state 
estimation protocols (like "minimal qubit tomography"), and the results are compared, discussed. 
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I. INTRODUCTION 

The main aim of this paper is to introduce a contin- 
uous positive operator valued measure (POVM) [11,3, 
which is equivariant under the symmetry group SO(3, M) 
of the qubit state space, and investigate the properties 
of the qubit state estimation protocol based on the gen- 
eralized measurement defined by this POVM. There are 
several similar qubit state tomography protocols in the 
literature [3, |3| operating with cleverly chosen discrete 
measurements, which respect only a restricted, discrete 
group of symmetries. Probably the most famous of them 
is the "minimal qubit tomography" protocol [5,]. I found 
it preferable to be a bit didactic here and to insert a brief 
account on these discrete qubit tomography protocols, for 
better comparability of the methods and results, and for 
the convenience of the reader. Thus the second section 
together with this introduction summarizes known re- 
sults, while the new investigations are put in section [1111 
The properties of the different protocols are compared in 
the concluding section. 

In quantum information theory Q the term quantum 
bit (or qubit) refers to the simplest nontrivial quantum- 
mechanical system, which has only two independent pure 
states. The Hilbert space of the qubit is C^, and its pos- 
sible states are described by the 2x2 density matrices 
p, which can be decomposed in terms of the three Pauli 
matrices 
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in the following way: 
1 
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Irl = r < 1. 



(The condition r < 1 ensures that p(r) > 0.) It means 
that the possible states of the qubit are in one-to-one 
correspondence with the points of the unit ball in M'^, 
which is conventionally called the Bloch ball. The surface 
of the ball, the Bloch sphere (r = 1) represents the pure 
states, i.e., rank one projectors. The mapping r i-^ p{r) 
respects convex combination (or "averaging"), so 

piar, + /3r2) - ap{r,) + ppir^); (p(r)) = p((r)) (3) 

provided that a + /3 = 1. 

The state space of the qubit has a nontrivial symme- 
try group determined by the automorphism group of its 
event lattice. In terms of density matrices the symmetry 
transformations are conjugations by unitaries, while in 
the Bloch ball picture the symmetries correspond to or- 
thogonal rotations of the Bloch ball. The unitary group 
SU{2, C) is a two-fold covering of the orthogonal group 
50(3,R), so they both have the same Lie-algebra. The 
structure of this Lie-algebra is reflected in the multipli- 
cation relations of the Pauli matrices, which we shall use 
in the sequel: 
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The principal problem of quantum state estimation, or 
quantum tomography is to give an accurate estimation 
p' for an unknown qubit state po by performing certain 
quantum measurements on multiple replicas of the un- 
known quantum bit. In the present approach it is essen- 
tial that all the replicas are in the same state pq, only one 
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measurement is performed on each single replica, and the 
choice of the measurement does not depend on the previ- 
ous results. Because of the inherent probabilistic nature 
of quantum mechanics the results of the measurements as 
well as the estimation p' itself are random variables. The 
estimation is unbiased if (p') = po, where (•) designates 
the expectation value taken over the possible outcomes 
of the measurements. The accuracy of the estimator p' 
is characterized by its variance ((P{p', {p'))), where d is 
an appropriate distance on the set of density matrices. 
A usual choice is the distance d^(pi,p2) = Tr(p2 — Pi)^ 
based on the Hilbert-Schmidt norm for selfadjoint matri- 
ces. 

It is worth noting that quantum state estimation is 
a much more delicate problem than its classical coun- 
terpart because of two main reasons: i) in quantum me- 
chanics two quantities usually cannot be measured at the 
same time with arbitrary high precision; and ii) measure- 
ments destroy the original state of the system. That is 
why in quantum tomography different, carefully chosen 
protocols exist for the measurements applied for the es- 
timation of the unknown state. 

In the second section we revisit three existing protocols 
and their basic properties. In all cases the measurements 
have a finite set of possible outcomes. The first protocol 
is based on von Neumann (projective) spin measurements 
in three orthogonal directions, thus there are three dif- 
ferent measurements, each having two possible outcomes. 
In the second protocol these three projective measure- 
ments are put together to obtain a (non-projective) mea- 
surement with six possible outcomes, based on a positive 
operator valued measure (POVM). The third one, the so 
called "minimal state tomography" protocol 5] is very 
much alike the previous protocol, but the number of pos- 
sible outcomes of the POVM measurement is reduced to 
four, which is a lower bound in qubit state tomography. 

In these three protocols the maximum likelihood 
method is applied to obtain the estimator p' . Unfortu- 
nately in certain cases the likelihood function may take 
its maximum value outside the Bloch-ball, what makes 
the further exact analysis rather complicated Dis- 
regarding this fact, we show that the "unrestricted" es- 
timator is unbiased (for the first and third protocol) or 
asymptotically unbiased (for the second protocol), and 
we calculate the variance of the (unrestricted) estimator 
in all cases. The variance goes to zero as the number TV 
of measurements is increased, what justifies that for large 
N and mixed states the unrestricted and exact estima- 
tors are essentially the same. For pure states, however, 
the unrestricted estimator is not a good choice. Another 
principal disadvantage of these discrete protocols is the 
fact that they do not respect the whole symmetry group 
of the qubit. 

This last defect is rectified in the third section, where 
we present a new protocol for qubit tomography based on 
a POVM, which is supported on the Bloch sphere, and 
equivariant under the symmetry group of the state space. 
An interesting novelty is that the POVM applied here is 



continuous^ i.e., the corresponding measurement has an 
infinite number of possible outcomes, namely all the pure 
states. Although the maximum likelihood estimator can- 
not be explicitly constructed, we present another simple 
unbiased estimator and calculate its variance. 

We conclude by comparing the results obtained for the 
different discrete and continuous qubit state estimation 
protocols. 

II. DISCRETE MEASUREMENTS 

Originally von Neumann defined quantum measure- 
ment as choosing one out of a set of pairwise complemen- 
tary events which form a complete system Translat- 
ing it into an algebraic language, a von Neumann type 
(or projective) measurement is defined by a complete set 
of orthogonal projections {Psjsgs. Here S is an appro- 
priate index set, and the projections satisfy the relations 

Y,Ps=I. Ps= P:, PsPr = S.^rPs- (5) 

ses 

In the state p the probability that the measurement re- 
sults in the event Pg is Ps = Tt{Psp). Usually it is con- 
venient to "label" the possible outcomes, i.e., the index 
set S by real numbers as and define a selfadjoint oper- 
ator A — X^ses '^sPs by its spectral decomposition. In 
this case we interpret A as an observable, the Os's are 
the possible values of the observable (measurement), and 
the expectation value of A is given by the well known for- 
mula Tr{Ap). But the essential part of the von Neumann 
type measurement is the orthogonal decomposition ([5]) of 
unity. Sometimes there is no natural way (or need) for 

the embedding S ^ R. In this case we can still speak 
about the probabilities ps, which give a classical proba- 
bility distribution on S, but (without further structure 
on S) the expectation value of the measurement has no 
meaning. 

This scheme can be generalized to the so called positive 

operator valued measure (POVM) p, g. In this case 

the identity is decomposed into the sum of (arbitrary) 
positive operators: 

J2Qs = I, Qs> 0. (6) 

The possible outcomes of the measurement are labelled 
by the index set S, and the probability of the outcome s 
in the state p is Ps = Tt(Qsp). As in the previous case, 
the Ps's define a classical probability distribution on S. 

The POVM measurements are also called weak mea- 
surements, particularly if the investigated system is cou- 
pled to another system and the POVM measurement is 
obtained from a projective measurement performed on 
the composite system [l|, Q . 

In this section we deal with discrete measurements, 
what means that the index set S is finite. The continu- 
ous case is conceptually the same, only technically more 
difficult. We turn to this in section Hill 
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Figure 1: The six projections belonging to the spin measure- 
ments in three orthogonal directions in the Bloch ball. 



A. Orthogonal spin measurements I. — Projective 

case 

In this protocol three different von Neumann spin mea- 
surements are performed in three orthogonal directions 
(see figure [TJ, so the projections are: 



P± = P(±Z) : 



2 



Py = p(±y) 



I±a. 



y 



(7) 



where x, y and z denote the unit vectors in the three 
orthogonal directions. 

Let N^^N+ + , Ny = N+ + Ny and = N+ + 
N~ be the number of spin measurements performed in 
the X, y and z direction, where is the number of 
results, and Nj^ is the number of results (fc = x, y, z). 
Furthermore, let po — p{ro) = ^{I+rocr) be the unknown 
true state of the qubit, and let p' = p(r') = i(/ + r'cr) 
denote its maximum likelihood estimator. 

The probabilities of the different measurement out- 
comes in the state po ^re 

p±^Tr(poP.±) = i±^, 

(8) 



i±yo 



l±zo 



so the likelihood function C{ro) for a given measurement 
statistics is 



which yields 



Ny 

N+-N- 



(11) 



A principal problem with this "unrestricted" estima- 
tor is that it may fall out of the Bloch ball. In this case 
either the true maximum place should be found on the 
boundary of the Bloch ball, or r' should be normed in 
some other way, e.g., by dividing it by |r'|. For the mo- 
ment we disregard this problem and continue as if r' were 
a true estimator. 

Using the binomial theorem we obtain that 
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in any of the three directions, and applying this and the 
previous equations ([8]) , ([9l) , pTjl a straightforward calcu- 
lation shows that the estimator r' is unbiased: 



r' = 



E 



JV+=0 Af+=OAf+=0 



r'/:(ro) = ro. 



(13) 



(The calculation can be carried out separately in the 
three spin directions.) 

Now let us determine the variance of the unrestricted 
estimator. Again a simple application of the binomial 
theorem shows that 

/iV+(iV 



/iV+(iV+ - 1)\ 
\ N{N-\) I 
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(14) 



in any of the three spin directions. The distance square 
between the expectation value (p') — po and the estima- 
tor is 



^(i-o) 
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(9) 



Varying the argument of the likelihood function, it takes 
the maximum value either on the boundary of the Bloch 
ball, or at the zero of its gradient: 



grad/:(r') =/:(r') 
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(10) 
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d'{p\ ip')) = Trip' {p')f ^ - Tr(Ar^)2 = -Ar^ 

(15) 

(Here Ar is a shorthand notation for the vector r' — (r'), 
and Ar is its absolute value.) Again the three directions 
decouple, what simplifies the calculations. The average 
of As^ can be determined with the help of (fT4| : 



1 



Nx 
(16) 

It means that in the states po = the variance in the 
X direction is zero, as it is expected. The total variance 
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is obtained by summing similar expressions for the a;, y 
and z direction: 



1 



1-^0 
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2iV, 



(17) 



If we assume that — Ny — — -j, where A'' = 
- Ny + is the total number of measurements, then 



the variance has the form 



3r§ 



2N 



(18) 



The variance tends to zero as the number of measure- 
ments goes to infinity, what justifies that for large N and 
for mixed states the unrestricted estimator pT|) practi- 
cally coincides with the true maximum likelihood estima- 
tor. 



B. Orthogonal spin measurements II. 



POVM 



This protocol is also based on the six projections P^y ^ 
appearing in the spectral decomposition of the three or- 
thogonal spin operators (see equation (O), but this time 
we form a single POVM out of them, by decomposing 
the unity in the following way: 

i = Qt + Q^+Qt + Qy+Qt + Q^, (19) 



where = ^P^ for k e {x,y,z} [see equation ([7])]. 
Performing this POVM measurement on each replica of 
the qubit, the six different results are obtained with the 



probabilities 



g± = Tr(poQ^) 



1 ± 3^0 

6 ' 



6 ' 



(20) 



provided that the qubit is in the state po = p{^a)- 

Let Nk — Nj^ + Nj^ (for k e {x, y, z}) be the number 
of measurements with any of the two possible results in 
the k direction, and let N — N^ + Ny + N^ denote the to- 
tal number of measurements. The probability of a given 
measurement statistics {N^} in the state po is described 
by the multinomial distribution: 



^(ro) 



m n 



N2 



(21) 



The unrestricted maximum likelihood estimator r' is 
obtained again by determining the zero of the gradient 
of £(r'), and we get formally the same expression pT|) as 
in the previous subsection. Now, however, an unavoid- 
able problem is the fact that N^, Ny or Nz may be zero 
with positive probability! If it happens then we have 
no information at all about certain component (s) of the 
Bloch vector r', and the division in (fTT|) is meaningless. 
In this case the most natural thing is to use the estima- 
tion for the appropriate component of r'. (We have 
to take care of this in the evaluation of the expectation 
values.) Of course this 'ad hoc' choice will bias the esti- 
mator towards zero. In addition, the estimator may fall 
out of the Bloch ball; we disregard it in the foUowings. 

Applying the multinomial theorem it is easy to show 
that 
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(1 - qx) 
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N 

Nx 



(1 - 9.)^-^^-<z^ 



(22) 



(Here qx = qt + Ix ~ \^ independently of the true state 
Po-) Similar expressions are valid in the other directions, 
and it follows that the estimator r' is biased, but asymp- 
totically unbiased: 

,23) 



Note that the factor 1— (2/3)^ is exactly the probability 
that none of the N measurement results lies in a partic- 
ular coordinate direction. 

To determine the variance of the estimator we need 
averages of the following type: 
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where 



71=1 

N 



n\2J 

^ AT / N ^ „ N 



(25) 



Here the last equation was obtained by integrating the 
fohowing sum with respect to q from to 



n=l 



dq ' — i \ n / n 



(1 + 9) 



N 



N-1 

E 

n=0 



(1 + g)". (26) 



Unfortunately the sum ((25|) defining cannot be fur- 
ther simplified, but it can be shown that lim7v-»cxD Fn — 
3. Indeed, the second sum in (|25|) diverges only logarith- 
mically, but it is multiplied with an exponentially de- 
scending factor, so it rapidly converges to zero. The first 
sum Sn — 

NJ2n=i ^2/3)^-" satisfies the recursion 



'JV-l-l 



S 



N 



N 



1, 



(27) 



which means that limjv^oo Sn = limTv^oo Fn — 3. 

Using and ([M]) the average of the error square of 
the X component of the estimator is 



Fjv 



2\N 



(1 -.?)+.?(!) i-(|) 



'2\^ 



(28) 



and similar formulas are valid in the other directions. By 
equation (jlSp the variance of the estimator is: 



Fjv (3 



,6 



N 



3rl 



2N 
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N 



2N 



(29) 



For large 1 this expression has the same asymp- 

totic behavior as the variance (fT5|) of the estimator using 
projective spin measurements in orthogonal directions. 



C. Minimal qubit tomography 

This protocol was introduced in The protocol is 
based on a POVM, which is minimal in the sense that 
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Figure 2: The four projections Pk at the vertexes at of a 
regular tetrahedron on the Bloch sphere. 

it contains only four positive operators {Qfej^^i in the 
decomposition of unity. These operators are constant 
multiples of the four projectors Pk = p{^k) being at the 
vertexes {afc}^^-^ of a regular tetrahedron on the Bloch 
sphere, as shown in figure [2j 
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1 
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(30) 



Qk ^p(afc) = 



(31) 



k=l 



Due to the tetrahedral symmetry, the vectors a^ have 
the following properties: 



1 4^ 
afea; = -- + -dkj 



(32) 



fe=i 



In the state po = p{ro) the probabihties pk of the four 
different measurement outcomes can easily be calculated 
with the help of the properties ^ of the Pauli matrices: 



Pk = Tr(poi9fc) 



Tr((/ + rocr)(/ + afccr)) _ 1 + rpafc 

8 ~ 4 ' 

(33) 

Denoting the number of fc-th measurement results by Nk 
and the total number of measurements by A'^ = X)fc=i ^k, 
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the probability of a given measurement statistics in the 
state pq is described again by the multinomial distribu- 
tion: 



>C(ro) 



iV! 



Ari!A^2!A^3!^4! 



P^'P^'-P^'P^^. (34) 



This likelihood function takes its maximal value at the 
zero of its gradient, which yields the equation 



1 + r'£ 

fc=i 



'afe 



(35) 



for the Bloch vector r' of the estimator p' = p(r'). 
For general at vectors it is hopeless to solve explicitly 
this equation, but for the symmetrically distributed vec- 
tors using their properties it is easy to check 
that the unique solution is 



k=l 



N 



(36) 



This solution is well defined for any measurement statis- 
tics {Nk}, but it may fall out of the Bloch ball. 

Using the multinomial theorem it is easy to show that 
{Nk) = Npk, so 

4 34 g . 4 . 

(r') = S^pfcafe - ^ afe + - f ^ afe o afc jro = ro, 

fc=i fe=i ^fc=i ^ 

(37) 

which means that the estimator p' is unbiased. 

Applying again the multinomial theorem, the formu- 
las (pS)) for the probabilities and the properties ([32| of 
the afc vectors, we obtain that 



4: 'i 



fc=i 



fc=i 
4 



^ {NiN - 1)pI + Npk 



(38) 



fc=i 
3 + r| 
12 



(iV^ „ TV) + iV. 



From this, using equations and (|32p we immediately 
get that 



fc=i 



9_rl 

N 



and finally, using (jlSp . for the variance we obtain: 



(r'2) - (r'>2 _ 9 - rg 



2iV 



(39) 



(40) 



Comparing it with the variances (flSl) and ((29)l obtained 
for the previous protocols, we see that in all cases the 
variances decrease as but for states close to pure 

states (i.e., if rg ^ 1) the minimal qubit tomography 
protocol is less performant by a factor of = |- 



III. CONTINUOUS, EQUIVARIANT 
MEASUREMENT 



Up to this point we have investigated three different 
qubit tomography protocols, which were based on (either 
projective or POVM) measurements with only a finite 
number of possible outcomes. In a finite d dimensional 
Hilbert space the maximal number of possible outcomes 
of a projective measurement is d, but this restriction does 
not hold for generalized, POVM measurements. In this 
section we construct a new qubit tomography protocol 
based on a continuous POVM measurement, which is 
equivariant under the symmetry group of the qubit sys- 
tem, and supported by the set of pure states, i.e., by the 
Bloch sphere. 

The attribute "continuous" refers to the fact that the 
investigated measurement has an infinite number of pos- 
sible outcomes, which form a topological manifold — the 
Bloch sphere in our case. The attribute "equivariant" 
refers to a kind of nice symmetry property, which de- 
serves a deeper investigation. 



A. Equi variance 

At an abstract level the symmetry group G of an (ei- 
ther classical or quantum) system is an automorphism 
group of its event lattice [1, 0]. States are probability 
measures on the event lattice, so G acts both on events 
and (by puUback construction) on states. In quantum 
mechanics events are projections in the Hilbert space 
Ti of the system, and by Gleason's theorem states 
are represented by density operators. By Wigner's the- 
orem 0, if dimTi > 2, every automorphism g of the 
projector lattice can be realized by conjugation with a 
unitary (or antiunitary) operator C/g, thus we obtain a 
left G-action on the lattice of projections (events) V, and 
a right G-action on the set of density operators (states) 
S: 



GxV 
SxG 



V. 



i9,P) 



gP = UgPU;; 

P9 UgpUg, 



(41) 
(42) 



and, by the puUback construction 



p{gP) - TiipUgPU;) = TiiU;pUgP) = {pg){P). (43) 

Performing many times a measurement with possible 
outcomes in the Borel space {S, B) we obtain a proba- 
bility distribution on the space S. Thus a measurement 
can be regarded as assigning to a quantum state p (i.e., 
probability function on V) a classical state Pp (i.e., a 
probability function on B). The simplest way to do this 
is pulling back p by a morphism P : B V between the 
classical and quantum event lattice: 



Pp = poP:B^ [0,1], pp{E) = Tr {pP{E)) . 



(44) 



7 



The classical quantum lattice morphism P is nothing 
but a projector valued measure, and the above formula is 
exactly the probability for the outcome of the projective 
measurement P being in E. Formally the generalized, 
POVM measurements are obtained by changing P to a 
positive operator valued measure Q on (S^B). 

Now let us assume that is a right G-space, and let 
7g : S" — > S" denote the transformation corresponding to 
g <E G. (So 7g o = Ihg ) Then its inverse image 
: B B defines a left action of G on We say that 
the POVM Q is equivariant with respect to the group 
actions 7"^ and U = U -U* (conjugation by U) if for any 
E € B and g € G the following diagram commutes: 



B 



B{n) 



B 



«1 

B{n) 



U,Q{E)U;=Q{^-\E)), (45) 



where B{'H) denotes the set of bounded operators in H. 

For projector valued measures this structure was dis- 
covered by George W. Mackey [13, [Hi, while he was 
studying the foundations of quantum mechanics [l3] , and 
he called it imprimitivity system. Mackey, in his im- 
primitivity theorem, classified and explicitly constructed 
the possible imprimitivity systems using induced repre- 
sentations. Since then induced representation became a 
very important tool in the theory of unitary group rep- 
resentation and harmonic analysis, while the imprimitiv- 
ity theorem turned out to be a cornerstone of quantum 
mechanics of free systems. Indeed, the equivariance con- 
dition ([IS]) , stated for projector valued measures on the 
three dimensional Euclidean space, is a very general and 
deep expression of the canonical commutation relations. 
Starting from this, many important properties (like the 
existence of mass, spin) of the elementary particles can be 
deduced, based on purely spacetime symmetry require- 
ments fTo|,[il[il,[i|. 

Borrowing the terminology from the projective case, 
we refer to the equivariance condition as general sys- 
tem of imprimitivity, more precisely, as the POVM Q is 
a generalized system of imprimitivity for U based on S. 
In the next subsection we give a very simple example for 
this in terms of a qubit. 



B. Equivariant POVM on the Bloch sphere 

First of all, we note that Wigner's theorem does not 
apply for two dimensional Hilbert spaces; the projector 
lattice of has a much bigger symmetry group than 
the symmetries induced by unitaries (and antiunitaries). 
People consider it as a pathological fact due to the low 
dimensionality of the space. Here we join this opinion 
and consider the group 5*0(3, R) — SU{2,C)/Z2 as the 
symmetry group G of the qubit. (The reason for the 
factorization is that conjugation with U and —U is the 
same transformation.) 



Furthermore, as the set 5* of possible measurement out- 
comes we choose the set of pure states, i.e., the Bloch 
sphere itself! It is a very clever choice for several reasons: 
i) By (H^ S* is a right G-space in a natural way. ii) Clas- 
sically pure states are Dirac measures on the phase space, 
thus in this respect S is the quantum analog of classical 
phase space. Is there any better measurement than the 
one returning the (or a?) phase space position (pure 
state) of the measured system? And finally Hi) pure 
states are rank one projections, i.e., positive operators, 
thus the POVM Q is "already there" on S, one has only 
to normalize it. 

The rank one projectors p{s) are parametrized by 
points s of the Bloch sphere (|s| = 1), unitaries have 
the form t/„ = e*"*^ (where |n| < 27r), and it can be 
shown that 



>(s)e- 



p(0-2„s), (46) 



where Om is the orthogonal rotation around the axis m 
at an angle |m|. (To derive this the equation [ncr, mcr] = 
2i(n X m)(T is needed, which is a straightforward conse- 
quence of (|3]) .) Thus in the Bloch sphere picture the sym- 
metry transformations are the rotations of the sphere! 
(We remark here that the nonunitary symmetries are ar- 
bitrary homeomorphisms of the Bloch sphere which map 
antipodal points to antipodal points.) Equations (|42p 
and (|46p also show that in this case the right G-action 
7n on the space S is 

7n(s) — representing point of U*p{s)Un = 02nS. (47) 

Now let CO be the normalized area measure on the 
unit sphere S, i.e., in spherical polar coordinates ("djip), 
dLj{d, If) = slRldddip, and for a Borel set E (Z S let the 
POVM Q be defined by 

Q{E) = 2 [ p(s)du;(s) = 2iu{E)p{J^^^ sd^(s)) , 
JseE 

(48a) 



or shortly 

dQ{s) = 2p(s)rfw(s), 



Q ~ 2puj. 



(48b) 



(The last equality in (|48ap is a consequence of ([3]).) It is 
clear that Q is normahzed, i.e., Q{S) — 2p{0) = I. 

Using the previous three formulas (|^5H15|) and the in- 
variance of the area measure uj = uj o under any 
rotation O, it is simple to show that Q satisfies the equiv- 
ariance condition (|^5)) : 

U^Q{E)U*^ = 2c^(^)p(/^^^0_2„sdc^(s)) 

= 2u:{E)p{J^^^_^^^sdw{s)) = Q(7n'(^))- (49) 

(It is worth noticing that here essentially the equa- 
tion and the invariance of lo was used.) 

We give the concrete formula of the POVM Q in spher- 
ical polar coordinates {'&,(p), although we will keep on 
using the abstract form (|48)) : 



dg(^,(p) 



sinz? 
47r 



1 -|-cosi9 e-'fsm-d 
e**^ sin I? 1 — COS'S 



dMip. (50) 
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Intuitively this POVM can be regarded as the hmit 
of discrete POVM's supported by more and more points 
scattered uniformly on the Bloch sphere. In the previous 
section we have seen examples for four (subsection III Cp 
and six fsubsection lll B[) supporting points with tetrahe- 
dral and hexagonal symmetry, respectively. 

Now let us investigate the distribution of the measure- 
ment Q provided that the system is in state po — p{^q)- 
The probability density function /o (with respect to lu) 
at the outcome s € 5 is 



Tr{poQ{dE)) 
uj{dE) 
2TV(p(ro)p(s)) = 1 



(51) 



i"oS, 



where dE is an infinitcsimally small area on the sphere 
S around s. (The last equality is a direct consequence 
of Q.) This probability distribution is unsharp, even for 
pure states! If rp = 1, then the maximum is /o(ro) = 2, 
and the minimum is /o(— Tq) = 0. 

Here it is worth noting that the POVM Q defined 
in (j48[) is not the unique solution of the equivariance 
condition (|45|l (but probably the most reasonable one). 
Indeed, for any a G [0, 1] the POVM 



Qo.{E) = aQ{E) + {1 - a)Lu{E)I 



(52) 



satisfies the equivariance condition (|45p. but the density 
function of Tr {poQa{-)) at s e 5 is 1 + arps, which is 
even less sharp then ([5T|) . 

Finally it is worth calculating the measure of a semi- 
sphere. The "center of mass" of the semisphere 5+(r) 
with midpoint r (|r| — 1) is at |, so by (|48ap 



Q((5+(r))=p(^ 



/ rcr 

2+- 



(53) 



Comparing it to the orthogonal projections (O, we see 
that if we are interested in a particular spin component 
of the state po, then it is better to perform a projective 
spin measurement than measuring Q. But what if we are 
interested in all components of the vector rg? 



C. Maximal qubit tomography 



this measurement statistics in the state po = p(vq) is 



N 



'C(ro) = n(l + i"onfe) 



(54) 



fe=i 



Its gradient at r' is C{r') i+l^n ; which yields the 

likelihood equation 



N 



^ l + r'i 



k=l 



rife 



(55) 



The second derivative of the logarithm of the likelihood 
function 



gy.2 



N 



nfe o rife 



(1 + rn,y 



(56) 



is everywhere negative definite, so the likelihood equa- 
tion (|55p has a unique solution for r', where C{r') is max- 
imal. Unfortunately for a general measurement statistics 
{^k}k=i ^^^^ solution cannot be analytically determined. 
(In subsection III CI a similar equation ([55)1 was obtained, 
which could be explicitly solved ([55]) because of the sym- 
metry of the fixed {afe}^^-^ vectors.) 

Instead of using the maximum likelihood estimator we 
introduce another obvious estimator 



N 



k=l 



with 



fiN) 



3 

N 



(57) 



where the (yet) unknown coefficient f{N) is determined 
from the expectational value of r'. For this aim we need 
two simple integrals: 



f nduj{n) — 0, f no ndLj{n) = —I 
JnGS Jnes 3 



(58) 



where lo is the normalized area on the unit sphere S. (The 
first integral is zero by symmetry, and the second integral 
can easily be calculated in spherical polar coordinates 
{^,ip), where dw(i?,^) = ^dMip.) 

Thus, using ([54|) . ([ST]), and then ([58|) . the expectational 
value of the estimator is 



In this subsection we investigate a quantum bit state 
estimation protocol based on the POVM Q : B{S) 
B{C'^) introduced in the previous subsection [see equa- 
tion We find it convenient to call this protocol 
maximal qubit tomography, since the whole set S of pure 
states constitute the possible measurement outcomes. 

Assume that performing N independent measurements 
on replicas of the qubit in the same state po, the measure- 
ment outcomes {nk}k=i C S are obtained, where each 
nfe represents a pure state on the Bloch sphere 5, i.e., 
nfe G M^, I nfe I = 1. Using ((5T|) . the likelihood function, 
i.e., the probability density function on of obtaining 



(r')=/(iV)V/ £(ro)nfed^c.(nfe) 

= Nf{N) f n(l + ron)dc^(n) = ^^^^ 
Jnes -J 



(59) 



which means that indeed, f{N) = should be chosen 
to get an unbiased estimation. (The integrals decouple, 
and for every k only the integral over nfe gives a nontrivial 
result, the other integrals are 1. The obtained result is 
also in accordance with 
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In order to calculate the variance we need another in- 
tegral, which can be easily obtained from ([55)1 : 

J J nm(l + ron)(l + rom)dtj(m)dw(n) = (60) 

neS meS 

Using this, we get that 

_ 9 + (A^ - l)rg 
N 

so, by (fT5|) . the variance of the (unrestricted) estima- 
tor dnZl) is 

2 (^^^ 

which exactly coincides with the one (|40p obtained for 
the minimal qubit tomography protocol. 

IV. CONCLUSIONS 

In table U we summarize the results obtained for the 
four different qubit state estimation protocols. (The vari- 
ance of the state p(r') is calculated from the Hilbcrt- 
Schmidt distance, and by equation ([TS]) it is half of the 
variance of r'.) 

It is clearly visible that all estimation schemes perform 
more or less equally well, although there are projective 
(*), discrete POVM (®, (g)), and continuous POVM mea- 
surements (o) among them. For states close to pure 
states (i.e., for ro ^ 1) the variance of the first two pro- 
tocols with hexagonal symmetry (* and ®) is a bit less 
than the variance of ® and O- On the other hand, the 
variance of the minimal and maximal qubit tomography 
protocol is exactly the same for all states. 

Based upon these examples we may draw some general 
conclusions on the unusual properties of POVM measure- 
ments. 

In classical Hamiltonian mechanics the phase space can 
be identified with the set of pure states of the system, 
which are simply Dirac measures concentrated on a single 
point of the phase space. This observation motivates to 
regard the set of pure states even in a finite dimensional 
quantum mechanical system as the quantum analog of 
phase space. 



A widespread paradigm of quantum mechanics, based 
on Heisenberg's uncertainty principle, is the fact that all 
the classical phase space variables cannot be measured 
simultaneously with arbitrary high precision. Originally 
this was stated for 'projective measurements, which are 
sharp in the sense that for every possible measurement 
outcome there is a state for which the specified outcome 
occurs with probability one. (Right after the measure- 
ment the system "jumps" into a state of this kind.) Co- 
ordinate and momentum operators do not commute, thus 
they do not possess a common projector valued measure. 

On the other hand, there are generalized POVM mea- 
surements which yield results corresponding to simulta- 
neous values of noncommuting operators! This fact is, 
however, not in contradiction with Heisenberg's uncer- 
tainty principle, since the result of the POVM measure- 
ment is unsharp, i.e., for any pure state of the system 
the measurement results have dispersed probability dis- 
tributions. A very simple example is the POVM mea- 
surement constructed in subsection IIII Bl which clearly 
demonstrates all the unusual features of generalized mea- 
surements. In summary, in contrast to projective mea- 
surements, a POVM measurement 

^ may have a continuous set of possible outcomes, 
even in a finite dimensional Hilbert space; 

may yield simultaneously values for incompatible 
(noncommuting) observables; 

but these simultaneous results are unsharp. (In 
particular, it means that repeated measurements 
do not give the same outcome.) 

— Furthermore, the information gains (^) and infor- 
mation losses (^) somehow compensate each other, 
so cleverly chosen POVM measurements yield more 
or less the same amount of information about the 
true state as projective measurements. 

^ Last but not least, POVM measurements can be 
much better adjusted to respect continuous symme- 
tries of the system than projective measurements. 

Finally we remark that POVM measurements (or weak 
measurements) are not pure mathematical constructions, 
they can also be experimentally realized by making an 
appropriate projective measurement on a composite sys- 
tem. 
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Protocol 


mean 


variance of p' 


* 


projective spin measurements in three orthogonal directions 


unbiased 


9 - 3rg 
2N 


® 


POVM measurement in six orthogonal directions 


asymptotically 
unbiased 


9 — 3rn 
2iV '''' 


minimal qubit tomography 
(POVM measurement in four tetrahedral directions) 


unbiased 


9-ri 
2N 


O 


maximal qubit tomography 
(POVM measurement uniformly on the whole Bloch sphere) 


unbiased 


9-rg 
2N 



Table I: The results obtained for the different qubit tomography protocols. 
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